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In pump-probe time and angle-resolved photoemission spectroscopy (TR-ARPES) experiments 
the presence of the pump pulse adds a new level of complexity to the photoemission process in 
comparison to conventional ARPES. This is evidenced by pump-induced vacuum space-charge effects 
and surface photovoltages, as well as multiple pump excitations due to internal reflections in the 
sample-substrate system. These processes can severely affect a correct interpretation of the data 
by masking the out-of-equilibrium electron dynamics intrinsic to the sample. In this study, we 
show that such effects indeed influence TR-ARPES data of graphene on a silicon carbide (SiC) 
substrate. In particular, we find a time- and laser fluence-dependent spectral shift and broadening 
of the acquired spectra, and unambiguously show the presence of a double pump excitation. The 
dynamics of these effects is slower than the electron dynamics in the graphene sample, thereby 
permitting us to deconvolve the signals in the time domain. Our results demonstrate that complex 
pump-related processes should always be considered in the experimental setup and data analysis. 


INTRODUCTION 

Angle-resolved photoemission spectroscopy (ARPES) 
offers a unique capability to investigate the electronic 
properties of quasiparticles in complex systems with en¬ 
ergy and momentum resolution. Recent developments 
have seen significant advancements in adding femtosec¬ 
ond time resolution to ARPES in a pump-probe scheme 
mM- New insights into the nature and origin of intrigu¬ 
ing physical phenomena, such as collective excitations 
and phase transitions, can be gained by observing the 
time-dependent changes directly in the band structure. 
More specifically, time-resolved ARPES (TR-ARPES) 
gives access to the time-dependent changes of the spec¬ 
tral weight in selected states in the Brillouin zone. The 
different time scales observed connect to the different en¬ 
ergy scales characteristic of the various kinds of interac¬ 
tions simultaneously present in the material [4]. Conse¬ 
quently, interactions felt by the quasiparticles in specific 
states can be disentangled by their characteristic time 
scales. This time, energy and momentum resolving ca¬ 
pability is a virtue precluded to time resolved all-optical 
techniques. It has opened unprecedented possibilities to 
study the quasiparticle dynamics in a wide range of con¬ 
densed matter systems, including the ultrafast amplitude 
oscillations in charge density wave materials [IHS], the 
quasiparticle recombination rate near the d-wave node 
in cuprate superconductors [DEHH], the dynamics taking 


place on the Dirac cone in graphene [SHISl and on the 
surface of topological insulators pAlII^ . 

A TR-ARPES experiment is performed by first excit¬ 
ing the system out of equilibrium by electron-hole pair 
creation using a femtosecond pump pulse with an energy 
below the sample work function. After a variable time 
delay, the excited electron distribution is photoemitted 
using a second probe pulse with an energy above the work 
function threshold, as illustrated in Fig. for the case of 
monolayer graphene. As the time delay is varied, the ac¬ 
quired ARPES spectra monitor the relaxation of this ex¬ 
cited distribution back to the equilibrium state. This re¬ 
laxation process involves a complex interplay of ultrafast 
electron-electron and electron-phonon scattering events. 
For monolayer graphene on SiC, several recent studies 
report that the initial thermalization of the excited car¬ 
riers to a hot Fermi-Dirac (ED) distribution can be ac¬ 
companied by a multiplication of the number of hot elec¬ 
trons with energies larger than the Fermi energy [12] [2QF 
EH, and that the subsequent cooling of the hot electron 
distribution happens on different time scales determined 
by energy relaxation via optical phonons and acoustic 
phonons, primarily assisted by disorder-scattering (so- 
called super collision processes) [9l [26ll28] . These pro¬ 
cesses are sketched in Fig. on the Dirac cones in movie 
snapshots. 

To extract this intrinsic electron dynamics from the ac¬ 
quired TR-ARPES data, great care should be exercised 
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FIG. 1: Schematic of a TR-ARPES experiment on graphene. Electrons (blue spheres) excited by an infrared pump pulse are 
photoemitted by an ultraviolet probe on an ultrashort timescale and detected by an electron analyzer. This provides energy-, 
momentum- and time-resolved snapshots of the electron dynamics around the Dirac cone. Such measurements provide direct 
information about electron and hole (yellow spheres) recombination via electron-electron (curled arrows in left Dirac cone 
sketch) and electron-phonon (wiggled arrows) processes. Interactions between electrons and optical phonons (middle Dirac 
cone sketch), acoustic phonons and impurities (right Dirac cone sketch) can be disentangled by the timescales revealed in the 
non-equilibrium signal. 


as several experimental effects might influence the mo¬ 
mentum and energy distributions of the photoelectrons 
before, during and after photoemission. In fact, in a low 
repetition rate laser source (< 1 kHz) the intensity of 
both pump and probe pulses is necessarily high in order 
to detect a pump-probe signal out of the noise within 
a reasonable experimental time window. This implies a 
high electron density in the cloud of photo-emitted elec¬ 
trons propagating towards the electron analyzer. The 
mutual Coulomb interactions between the photoelectrons 
create a vacuum space-charge effect that may lead to a 
shift of the electrons’ kinetic energies and a severe broad¬ 
ening in energy that can easily exceed the resolution set 
by the light source and the analyzer [29H3T] . Moreover, 
a very recent time-resolved x-ray photoelectron spec¬ 
troscopy study by L.P. Oloff et al. [32] shows that pho¬ 
toelectrons emitted by the x-ray probe beam may also 
be influenced by the electric field originating from sec¬ 
ondary electrons generated by the infrared pump. This 
may appear counterintuitive as the pump energy is typ¬ 
ically lower than the sample work function (< 4 eV). 
However, high-order nonlinear photoemission processes, 
where an electron simultaneously absorbs several pho¬ 
tons, can provide sufficient energy for the electrons to es¬ 
cape into vacuum [32l [33] . These processes can dominate 
as the pump power (^ 1 mW) is typically several orders 
of magnitude higher than the probe power (^ 1 nW). 
For semiconductors, an additional effect, originating from 
a light-induced voltage difference in the surface region, 
can affect the propagating photoelectrons. Specifically, 
at the surface, the photo-excited electron-hole pairs can 


be spatially separated by the band bending. This sepa¬ 
ration generates a transient electric field that affects the 
photoelectron outside the sample. Despite having been 
known for several years, it was only very recently that 
the implications of this surface photovoltage (SPV) for 
TR-ARPES experiments were addressed theoretically by 
S. Tanaka [34] and experimentally by S.-L. Yang et al 
on semiconducting GaAs [35] . 

In this paper, we investigate the extent of these per¬ 
turbing effects in the acquisition of TR-ARPES data 
from monolayer graphene on a semiconducting SiC sub¬ 
strate. We find strong indications of pump-induced elec¬ 
tric fields and pump beam reflections at the backside of 
the SiC substrate that affect the measurement on picosec¬ 
ond timescales. A thorough understanding of these ef¬ 
fects is important in order to suppress their impact on 
the experiment or deconvolve them from the intrinsic 
electron dynamics in the data analysis. 

EXPERIMENT 

The graphene sample investigated in this work was 
produced following a well-documented synthesis method 
that provides high-quality monolayer graphene on the Si- 
terminated face of SiC(0001). The graphene layer was 
decoupled from the substrate by hydrogen intercalation 
[Ml ET] such that the structural and electronic prop¬ 
erties probed by transport, Raman and ARPES mea¬ 
surements closely approached pristine graphene [37f[4Q] . 
The graphene is p-doped with a carrier concentration of 
« 5 X 10^^ cm“^ that places the Dirac point 240 meV 
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above the Fermi energy. The doping mechanism is caused 
by a spontaneous polarization in the bulk of the semicon¬ 
ducting SiC substrate m- The sample was cleaned by 
annealing to 500 K in ultra-high vacuum in order to re¬ 
move adsorbed water, and was held at room temperature 
throughout the entire experiment. The thickness of the 
SiC wafer was 390 fim. 

All TR-ARPES spectra presented in this work were ac¬ 
quired at the Artemis facility, Rutherford Appleton Lab¬ 
oratory using the TR-ARPES end-station da SSI. The 
extreme ultraviolet, ultrafast probe pulses provided at 
this facility by high-harmonic generation enable acquisi¬ 
tion of TR-ARPES data from the ^-point at the corner 
of the hexagonal Brillouin zone of graphene. The sample 
was excited by using a beam of 30 fs light pulses with a 
repetition rate of 1 kHz, a photon energy of hz/ = 1.55 eV 
and with fluences ranging from 1 mJ/cm^ to 6.6 mJ/cm^. 
The band structure around the Dirac cone was measured 
by a time delayed probe pulse of = 21 eV correspond¬ 
ing to the 13th harmonic of the laser fundamental. Eor 
the measurements in this work, a large time delay win¬ 
dow spanning from —100 ps to 100 ps was used in all 
scans. Negative time delays correspond to a situation 
where the probe pulse arrives before the pump pulse. 
The grid of time delay points was varied in such a way 
that a coarse grid was used for long time delays, while a 
fine grid was used to properly capture the hot electron 
signal around the time of the optical excitation that we 
define as time zero (t = 0). The time and angular resolu¬ 
tions were better than 40 fs and 0.3 °, respectively. The 
intrinsic energy resolution was 380 meV in these experi¬ 
ments. Additional fluence dependent broadening effects 
were present, leading to an apparent energy resolution 
denoted as AE, which we discuss in the following sec¬ 
tions. 


RESULTS AND DISCUSSION 
Snapshots of electron dynamics 

TR-ARPES measurements of the low-energy Dirac 
spectrum of graphene at a pump fluence of 2.2 mJ/cm^ 
are shown in Eig. The acquired dynamics is di¬ 

vided into four temporal regions labeled I, II, III and 
IV that are marked in Eig. |^a). Snapshots of the Dirac 
cone of graphene at selected time delays are presented 
in Fig. I^b), and corresponding intensity difference plots 
are shown in Eig. |^c). The latter are arbitrarily de¬ 
termined by averaging over all spectra taken in the time 
interval from —1000 fs to —100 fs and subtracting this 
average from the spectrum at any given time delay. The 
Eermi level is referenced to the spectrum at -100 ps. The 
intensity difference in Eig. [^a) was integrated over the 
boxed region in the intensity difference plots shown in 
Fig-i c). Time zero is fixed on the middle of the rising 


edge of the initial optical excitation (in region II) of the 
electrons in the sample. 

Several unexpected features emerge in the dynamics of 
the intensity difference in Eig. [^a). Perhaps most strik¬ 
ing is the presence of two large peaks (one in region II and 
one in region III) both originating from a pump excita¬ 
tion. The first occurs as the pump beam impinges on the 
surface and the second is due to a reflection of the pump 
beam at the back of the substrate. We present a more 
detailed discussion of this below. Eurthermore, in the 
temporal regions I and IV a depletion of intensity differ¬ 
ence with respect to the averaged spectrum immediately 
before the excitation is seen towards both large negative 
and large positive time delays. Comparing the spectrum 
at a time delay of —100 ps to the one at —540 fs (here 
defined as the boundary between regions I and II) in Eig. 
Hb), one observes that the spectrum at -100 ps is shifted 
rigidly towards higher binding energies and appears less 
broad. The intensity difference in Eig. [^c) at -100 ps 
reflects this behavior directly as there is a surplus (deple¬ 
tion) of intensity towards higher (lower) binding energies. 
In region IV it is found that the spectrum also shifts to 
higher binding energies with a small shift already evident 
around 10 ps. We can therefore define an apparent chem¬ 
ical potential fi and apparent energy resolution AE that 
both change with time. 

Optical excitation of electrons is observed in region II 
as a large peak, which rapidly decays on a femtosecond 
time scale and is fully relaxed after 6 ps. This behavior 
reflects hot electron dynamics in graphene i.e. thermal 
relaxation of hot electrons by phonon scattering. The 
additional pump excitation triggers this dynamics again 
in region III. The fact that the electronic system is heated 
twice is also evident from the smearing of the ED edge in 
the spectra and difference plots at both 40 fs and 6.84 ps 
in Figs. [2];b)-(c). 

Influence of pump-induced electric fields 

In order to track more closely the evolution of /i and 
the change of AE, we apply a procedure that accurately 
provides the energy distribution curves of the electrons 
at all time delays by analyzing momentum distribution 
curves (MDCs) of each spectrum as described in detail 
elsewhere [44]. Examples of such MDC-derived energy 
distribution curves and the fits to a ED function convo¬ 
luted with a Gaussian energy function with width AE 
are shown for several time delays in Eigs. ia)-(b). 

The complete time dependence of /i in regions I and IV 
is given in Eig. [^c). In region I the value of /i moves to 
lower binding energies until a plateau is reached around 
the time of the optical excitation. In region IV fi moves 
back to higher binding energies. A similar trend is ob¬ 
served for AE in Eig. I^f). In regions II and III the 
position of /i changes due to the hot electron dynamics 
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FIG. 2: Electron dynamics in graphene induced by a pump laser pulse with an energy of 1.55 eV and a fluence of 2.2 mJ/cm^ 
over a ±100 ps time window: (a) Time dependent intensity difference summed within the boxed region given in the intensity 
difference plots in (c). The dynamics has been divided into regions I, II, III and IV with representative TR-ARPES spectra 
shown in (b) for each region and for the time delays defining the boundaries between the regions, (c) Intensity difference 
obtained by subtracting a spectrum averaged over the time interval from —1000 fs to —100 fs from the spectrum at the given 
time delay. 


intrinsic to the sample. The population of photoinduced 
electrons above the Dirac point forces the chemical po¬ 
tential towards lower binding energies as seen in the fits 
in Figs. I^d)-(e). This behavior is perfectly consistent 
with a simple calculation of /i for graphene with elevated 
electronic temperatures assuming that the charge density 
is conserved [12]. 

The behavior of fi and /S.E in regions I and IV does not 
reflect the instrinsic electron dynamics of the graphene 
but is the result of pump-induced external electric fields 
that affect the kinetic energy of each photoelectron emit¬ 
ted by the XUV probe pulse as it propagates to the elec¬ 
tron analyzer [294I3T] . As mentioned above, such fields 
may originate from either a vacuum space-charge effect, 
a SPV effect in the SiC substrate |35| or a combina¬ 
tion of both. The space-charge is caused by electrostatic 
interactions between the probe photoelectrons and sec¬ 
ondary electrons excited by the pump (1.55 eV) to ener¬ 


gies higher than the sample work function by nonlinear 
multiple photon absorption processes. These perturbing 
effects, which can persist on the nanosecond timescale, 
are particularly effective at the low repetition rate and 
high fluences we use, which yield large energies per pulse 
j32l l33] . This type of nonlinear process can also lead to 
the creation of electron-hole pairs in SiC by exciting va¬ 
lence band electrons across the band gap. A SPV is then 
created, generating a dipole field that will also affect the 
kinetic energy of the probe photoelectrons jSH |35] . The 
sketch in Fig. |^i) illustrates the event sequence for these 
effects at positive delay times (t > 0) corresponding to 
the dynamics in region IV. Initially the pump simultane¬ 
ously induces a SPV and generates secondary electrons 
that propagate from the sample surface with low kinetic 
energies. When the probe arrives at a later time de¬ 
lay, it photoemits electrons which are decelerated as they 
move towards the secondary electron charge cloud. Due 
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FIG. 3: (a)-(b) MDC-derived energy distribution curves (dots) at the given time delays obtained from MDC fits to the spectra 
and htted to a FD function convoluted with a Gaussian (lines). The stated values for the apparent chemical potential and 
the apparent energy resolution /\E are results of the hts. (c)-(e) Fitted position of fi for (c) long positive (region IV) and 
negative (region I) delay times and for (d)-(e) time-domain regions II and III where the pump-induced hot electron dynamics 
in graphene takes place, (f) Fitted values of /\E as a function of time delay. The energy broadening is not dehned in regions II 
and III where the electronic system is heated. The pump laser fluence in (a)-(f) is 2.2 mJ/cm^. (g) Apparent energy resolution 
at a time delay of -200 fs and (h) chemical potential shift A/i defined in (c) as a function of pump laser fluence. (i)-(j) Sketch 
of pump-probe scenarios at (i) long positive delays (region IV) and (j) long negative delays (region I). Blue spheres signify 
photoelectrons emitted by the probe and red spheres signify secondary electrons generated by the pump. The red disk is a 
surface photovoltage (SPV) induced in the SiG substrate by the pump pulse. The time it takes a photoelectron to be detected 
at the analyzer is much longer than the typical time delay between the pump and probe pulses. 


to the much higher kinetic energies of the photoelectrons 
emitted by the UV pulse, they pass through the charge 
cloud and are re-accelerated towards the electron ana¬ 
lyzer. Photoelectrons emitted by the probe at a negative 
time delay (t < 0) in region I of our data will also be 
affected by these pump-induced electric fields since these 
will be generated while the photoelectrons propagate to¬ 
wards the analyzer as sketched in Fig. |^j). In both cases 
the acceleration of the probe photoelectrons as they prop¬ 
agate in vacuum causes a rigid shift in kinetic energy as 
well as a broadening in their energy distribution, as ob¬ 
served in Figs, [^c) and 3(f), respectively. We note here 
that a distinction between SPV and vacuum space-charge 


is difficult from our data. The similar slow dynamics in 
both regions I and IV suggest that space-charge is the 
dominating effect since the SPV should diminish within 
a few picoseconds at positive delay times due to fast re¬ 
combination of electron-hole pairs in the semiconducting 
substrate [35]. Furthermore, SiC has an indirect band 
gap of 3 eV such that electron-hole pair separation 
by the 1.55 eV pump pulse would require both multi¬ 
ple photon absorption and transfer of momentum by e.g. 
phonons, which seems unlikely. In the experiments on 
the observation of the SPV effect in GaAs by S.-L. Yang 
et a/., the effect of space-charge could be ruled out due to 
the very high repetition rate of the laser (80 MHz) and 
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a low pump fluence (25 /iJ/cm^) [35] in contrast to our 
experiments. 

The change of energy broadening with time, shown in 
Fig. if), leads to an estimate of the inherent energy 
resolution of the experiment of 380 meV. This value is 
found at long time delays before the energy distributions 
of the photoemitted electrons are severely affected by the 
abovementioned pump-induced electric fields. The ap¬ 
parent energy resolution during the intrinsic hot electron 
signal in graphene corresponds to the plateau reached 
in Fig. if) just before the optical excitation, where 
the pump-induced energy broadening is most severe. 
Throughout the dynamics in regions II and III it is as¬ 
sumed that this energy broadening remains fixed, and 
changes in the width of the Fermi edge are fully deter¬ 
mined by the variation of the electronic temperature T^. 

The observed energy shifts and broadenings are a func¬ 
tion of applied pump fluence as shown in Figs. Hg)-(h). 
Here, the value for /S.E is determined from the plateau 
just before the arrival of the optical excitation. Below 
1 mJ/cm^ the apparent energy resolution is similar to 
the intrinsic energy resolution we could determine in this 
experiment, but above this fluence broadening rapidly 
sets in. Above 4 mJ/cm^ the broadening is so large that 
it becomes difficult to discern spectroscopic features. The 
related chemical potential shift A/i defined in Fig. [^c) 
shows that the long time delay dynamics in regions I and 
IV becomes negligible at low fluences, while extensive 
shifts of up to 1 eV occur above 6 mJ/cm^. This be¬ 
havior is consistent with the fact that the multiple pho¬ 
ton absorption processes leading to either space-charge or 
SPV become more efficient with higher fluence. This flu¬ 
ence dependence is highly nonlinear as seen in the data 
in Figs. I^g)-(h). We also observe in the experiments 
that the energy shift and broadening are very sensitive 
to both beam spot size and the illuminated area on the 
sample, which relate to the sharpness of the beam profile 
and the sample morphology. It is therefore crucial in any 
experiment to keep the fluence sufficiently low that these 
effects do not overshadow the intrinsc dynamics in the 
sample. 

Double pump excitation of graphene 

The electron dynamics in regions II and III reflect the 
hot electron dynamics in graphene discussed in great de¬ 
tail elsewhere mm- The surprise here is that it occurs 
twice and with nearly the same intensity. The MDC- 
derived energy distribution curves in Figs. |^a)-(b) have 
been fitted by FD functions and a Gaussian resolution 
function with AE set by the apparent resolution imme¬ 
diately before the arrival of the first optical excitation. 
Since /i and Tg are free parameters, it is possible to un¬ 
ambiguously extract the time dependence of Tg for the 
hot electrons in graphene. For each Tg trace in regions 


(a) (b) 



Gps) 

FIG. 4: (a)-(b) MDC-derived energy distribution curves for 
(a) the first optical excitation and (b) the second excitation. 
Markers correspond to data at the given time delays and lines 
are fits to FD functions with the given electronic temperature 
convoluted with a Gaussian resolution function, (c) Fitted 
electronic temperature for first (round markers) and second 
(diamond markers) excitations in regions II and III, respec¬ 
tively. The red lines are exponential function fits with the 
given time constants for the decay parts. The time difference 
between the onset of the two excitations is 6.8 ps as marked 
by the double headed arrow. 


II and III, displayed in Fig. [^c), the rising edge and the 
subsequent decay are fitted by a product of exponential 
functions. The decay part is fitted by two exponential 
decays, which provide the time constants stated in Fig. 
I^c). Maximum electronic temperatures of 3200 K and 
3000 K are reached with a time difference of 6.8 ps. In 
both cases a decay with two similar time constants sets in. 
The fast decay 250 fs) is caused by optical phonons, 
while the slow decay (~ 2000 fs) is caused by the joint 
effect of acoustic phonons and impurities ElIIH]. 

We are able to confirm that the occurrence of two peaks 
originates from a back reflection of the pump beam in 
the substrate by using the measured time delay between 
the two excitations. The first excitation occurs once the 
pump reaches the sample, while the second excitation 
takes place after the pump has traveled back and forth 
in the SiG crystal. Using a refractive index of 2.6 for 
SiC at 800 nm, the measured delay of 6.8 ps would cor¬ 
respond to a thickness of ~ 390 /im of the sample, which 
matches the actual thickness of the crystal. By using a 
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sample with a different thickness of the SiC substrate, we 
have checked that the time delay between the two excita¬ 
tions changed accordingly. The fact that the maximum 
electronic temperature is high for both excitations shows 
that very little energy of the pump is dissipated during 
its initial interaction with both graphene and SiC. Since 
two excitations occur, the space-charge and SPV effects 
are also induced twice. It is difficult to predict how this 
affects the overall dynamics, and this could play a role 
for the slightly different decay times extracted for the two 
excitations in Fig. [^c). 

SUMMARY AND PERSPECTIVES 

Summarizing, we have shown that TR-ARPES spectra 
of the Dirac cone in monolayer graphene on SiC exhibit 
strong shifts and broadening in energy. The observed 
temporal evolution of these perturbations are ascribed 
to a complex interplay of vacuum space-charge as well 
as a surface photovoltage in the SiC substrate, induced 
by the pump beam via multiple photon absorption pro¬ 
cesses. Additionally, the pump beam is found to be re¬ 
flected by the back surface of the substrate, giving rise 
to a duplicate optical excitation signal of hot electrons in 
graphene. 

With the present experimental setup for TR-ARPES 
such effects are most easily avoided by lowering the flu- 
ence of the pump laser into a regime where the SPV and 
space-charge effects are negligible compared to the in¬ 
trinsic electron dynamics of the sample. A further explo¬ 
ration of the SPV effect, e.g. to disentangle its contribu¬ 
tion to the dynamics on the long picosecond time scales 
from the photoelectron space-charge cloud, would require 
a substrate with a different band gap. Alternatively, one 
could use graphene on a metal substrate to avoid SPV al¬ 
together, but this reduces the pump-probe signal due to 
screening of the pump laser field by the metal [13] . These 
issues with pump (and probe) induced external electric 
fields that interfere with the photoemission process are 
expected to become less of a problem with future TR- 
ARPES setups utilizing laser sources with much higher 
repetition rates since the energy per pulse will be signif¬ 
icantly reduced. 
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